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\ 

^ This  report  describes  both  theoretical  and  experimental  research 
performed  on  the  subject  of  slender  wing-body  configurations  at  high 
angles  of  attack*,  under  Grant  No.  DAERO-78-G-021  during  the  period 
October  1977  - November  1978. 

-•The  theoretical  part  of  the  research  consists  of  a method  for 
the  computation^  of  longitudinal  aerodynamic  characteristics  of  com- 
plex wing  configurations  at  high  angles  of  attack,  and  in  subsonic  flow 
including  non-linear  vortex  lift  caused  by  leading-edge  or  side-edge 
separations . 

The  modified  vortex  method  developed  can  handle  closely  coupled 
canard  configurations,  dihedral  angles,  flaps,  "saw-tooth" ,strakes , 
elevators,  as  well  as,  ground  effect  problems.  An  outline  is  made  for 
including  a body  and  for  calculating  the  influence  of  vortices  shed  from 
bodies  at  high  angles  of  attack,  in  the  present  method. 

Some  of  the  calculated  results  are  presented,  and  those  compare 
favourably  with  experimental  data. 

•The  experimental  part  of  the  research  consists  of  a study  made 
on  the  effects  of  jets  injected  near  the  nose  of  a body,  on  the  side 
forces  and  yawing  moments  appearing  at  high  angles  of  attack  at  sub- 
sonic and  transonic  speeds . 

One  of  the  main  conclusions  from  the  results  of  the  experiments 
is  that  a great  reduction  of  side-force  and  yawing  moment  can  be 
achieved  using  small  amounts  of  air  injection  by  this  method.,.  Further 
study  is  need  especially  within  the  high  subsonic  and  the  transonic 
regions . 
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LIST  OF  SYMBOLS 
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aspect  ratio,  b /S 

wing  span 

chord 

induced  drag  coefficient  D/qS 

lift  coefficient,  L/qS 

section  lift  coefficient 

pitching  moment  coefficient  M/qS  C , 

r®  r 

yawing  moment  coefficient  N/qS  C 

blowing  coefficient  tn.u./i  p U2S 

i .1  2 « 

side  force 
drag  force 

geometric  influence  matrix 
index 

index  (or  subscript  for  "let") 
strength  of  a vortex  segment 
lift  force 

nozzle  mass  flow  rate 
pitching  moment,  or  Mach  number 
yawing  moment 

number  of  subdivisions,  chordwise 
number  of  subdivisions,  spanwlse 
dynamic  pressure,  7>V 

distance  from  a point  to  a vortex  segment 
cross  flow  Reynolds  number 
wing  area,  or  reference  area 
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theoretical  let  velocity  (assuming  fully  expanded  isentropic 


flow) 


freestream  velocity 


component  of  free  stream  velocity  normal  to  a surface 
velocity  disturbances  in  x,  y,  z directions  accordingly 


Cartesian  coordinates 


position  of  center  of  pressure 


angle  of  attack 


air  density 


geometrical  angles 
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This  report  summarizes  the  research  performed  under  Grant  No. 
DAERO-78-G-021  during  the  period  October  1st  1977  to  November  30th 
1978. 

High  angle  of  attack  performance  is  one  of  the  most  important 
regions  in  the  flight  envelope  of  the  new  generation  of  aircraft 
and  missiles.  The  flow  phenomena  Involved  in  high  angle  of  attack 
flight  is  very  complicated  and  therefore  a most  challenging  area 
for  the  aerodynamic  research  effort.  The  analysis  may  involve  com- 
plex configurations  of  wings  and  bodies  on  which  flow  separation 
occur.  In  such  flows  strong  vortices  may  be  generated  causing  non- 
linear behaviour  of  the  aerodynamic  coefficients.  Furthermore, 
configurations  with  strong  interaction  between  closely  coupled  lift- 
ing surfaces  and  their  wakes  may  be  analyzed.  Another  well  known 
phenomenon  is  the  appearance  of  relatively  large  side  force  and 
yawing  moments  acting  on  bodies  at  high  angles  of  attack. 

Some  of  the  phenomena  involve  strong  viscous  effects  such  as 
the  not  yet  fully  understood  "vortex  break-down"  phenomenon,  occur- 
ing  at  high  angles  of  attack  when  strong  vortices  are  involved. 

The  present  efforts  of  aerodynamicists  to  analyze  these  problems 
are  concentrated  on  developing  theoretical  models  and  computational 

me  thods . 

This  report  describes  the  application  of  the  computational  method 
developed  by  Rom,  Almosnino  and  Zorea  [18]  , [17]  for  obtaining 


the  aerodynamic  characteristics  of  complex  configurations  at  high 
angles  of  attack  to  slender  bodies  with  or  without  low  aspect  ratio 
wings . 

Some  experimental  work  carried  out  by  Rom,  Sharir  and  Portnoy 
[45]  studying  the  effects  of  lets  inlection  on  the  reduction  of  side 
forces  acting  on  a slender  body  at  angles  of  attack  is  presented  as 
indication  to  methods  of  control  of  the  aerodynamic  effects  at  these 
flipjht  conditions. 
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CALCULATION  OF  THE  AERODYNAMIC  CHARACTERISTICS  OF  SLENDER  WINGS 
AND  BODIES  AT  RICH  ANGLES  OF  ATTACK. 


There  ere  a number  of  methods  for  the  calculation  of  aerodynamic 
characteristics  of  wings  and  bodies,  presented  by  Maskew  [1,2],  Margason 
and  Lamar  [3],  Woodward  [4],  as  well  as  by  the  Boeing  method  [5,6],  How- 
ever, those  methods  are  limited  in  their  solutions.  For  exanple,  Maskew' s 
method  [1,2]  includes  the  effect  of  the  rolling  up  of  the  wake,  but  is  un- 
able to  model  the  leading  edge  vortices  and  the  non  linear  lift  generated 
due  to  their  effects.  The  model  of  the  body  in  that  method  is  unable  to 
represent  the  vortices  shed  from  the  nose  or  from  other  parts  of  the  body 
which  cause  non  linear  effects.  Margason  and  Lamar  [3]  or  Woodward's 
program  [4]  do  not  take  into  account  the  rolling  up  of  the  wake,  and 
therefore  cannot  evaluate  the  aerodynamic  characteristics  at  high  angles 
of  attack,  for  lew  aspect  ratio  wings  with  side  edge  or  leading  edge 
vortex  separation.  The  first  program  [3]  represents  the  body  in  a 
too  simplified  manner,  taking  instead  only  its  prolection  on  a hori- 
zontal plane,  while  the  latter  [4]  represents  the  body  by  source  panels, 
which  are  a fair  representation  for  thickness,  but  are  unable  to  model 
vortex  separation,  thus  failing  at  moderate  and  high  angles  of  attack. 
Modification  of  the  last  two  methods  try  to  evaluate  the  nonlinear  lift 
by  superimposing  Polhamus  [7]  and  modified  Polhamus  [8]  methods  on  the 
linear  results.  However,  in  these  cases  the  results  are  not  always 
accurate  and  do  not  model  the  flow  field,  the  local  effects,  the  cor- 
rect load  distribution  or  effects  of  close  Interactions. 

The  method  developed  by  the  Rubbert  [5,6]  et  al . using  a doublet 
lattice  (which  is  equivalent  to  the  vortex  lattice)  resulted  in  a much 
refined  program,  and  its  results  are  accuracte  in  many  cases.  The 
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method  takes  into  account  the  rolling  up  of  the  wake  and  its  effect  on 
the  wing.  However,  the  use  of  Smith's  [91  model  of  conical  flow  to 
describe  the  initial  form  of  the  leading  edge  vortices  results  in  an 
lnconq>lete  model  which  cannot  deal  with  secondary  separations,  and  which 
seems  too  complicated  for  the  description  of  close  interaction  between 
the  vortex  flow  field  and  wakes  of  a closely  coupled  canard-wing  con- 
figuration. The  description  of  the  body  in  this  method  lacks  the 
possibility  of  representing  body  shedding  of  vortices.  The  problems 
discussed  here  are  also  typical  to  other  methods  which  try  to  solve 
conq>lex  configurations  of  general  form  [10,  11,  12,  19,  20,  21,  22,  23, 


In  the  NASA  Workshop  (Ref.  1),  the  papers  dealing  with  configu- 
rations which  Include  bodies  use  the  ideal  of  "Annular  Wing"  formed 
by  a vortex  lattice,  (as  for  example  appearing  in  the  works  of  Thomas .Miranda, 
Palko,  Heltsley  in  Ref.  1).  However,  none  of  these  works  deal  with 
the  rolling  up  of  the  wake,  nor  with  shedding  vortices  from  the  body. 

The  shedding  of  a pair  of  vortices  from  slander  bodies  may  occur 
already  at  quite  small  angles  of  attack.  The  pair  of  vortices  remain 
symmetrical  up  to  about  25-35  degrees,  and  then  occurs  a phenomenon 
of  asymmetry  in  the  flow. 

Description  of  the  phenomena  of  body  vortices  may  be  found  for 
example  in  the  works  of  Mendehall  and  Nielsen  [25],  Clark  and  Nelson 
[26],  Angellucci  [19],  [20]  and  Smith  [46].  A review  of  the  method 
trying  to  describe  bodies  with  vortex  shedding  appears  in  a paper 
by  Nielsen  [27],  Further  Investigations  of  this  problem  are  discussed 
by  Angelucci  [19],  [20],  Fiddler  [28],  and  Marshall  and  Deffenbough  [29]. 


Most  of  these  works  use  semi-empirical  data  for  calculations. 

2 .1 , The  Non-linear  Modified  Vortex  Lattice  Method  Calculation  of 

Longitudinal  Aerodynamic  Coefficients  of  Slender  Wings  at  High 
Angles  of  Attack 

Rom  and  Zorea  [13,14,  15,  16]  Introduced  a modified  vortex  lattice 
method,  handling  the  wing  and  its  wake  as  one  system,  enabling  the 
shedding  of  vortices  not  only  from  the  trailing  edge  or  the  side-edge, 
but  also  from  the  leading  edge  or  from  within  the  wing  surface.  The 
vortices  which  were  divided  into  a number  of  finite  segments,  were 
allowed  to  roll  up  and  form  the  structure  of  the  three  dimensional 
wake.  Very  accurate  results  were  obtained  for  various  rectangular, 
triangular  and  swept-back  low  aspect  ratio  wings,  including  their 
non  linear  characteristics.  Special  care  was  taken  for  the  treat- 
ment of  the  singular  vortex  element  to  avoid  numerical  problems,  and 
to  represent  closely  the  actual  phenomenon. 

This  research  program  which  was  supported  in  part  by  AFOSR  Grant 
No.  71-2145  has  been  expanded  to  include  the  body  effects.  Some 
results  of  this  new  effort  is  being  supported  by  DAERO  Grant  No. 
78-G-C21.  This  further  development  is  presented  by  Rom,  Zorea  and 
Almosnlno  [17,18],  where  a general  method  for  calculation  of  the  aero- 
dynamic characteristics  of  wings  and  multiple  lifting  surfaces  in  sub- 
sonic flow  at  high  angles  of  attack  was  developed.  The  new  method 
can  handle  coi^>lex  pi  an  forme  such  as  closely  coupled  canard-^rlng-tail 
combinations,  various  flaps,  elevators,  winglets,  dihedral  angles  as 
well  as  ground  effects.  This  presentation  of  the  analysis  is  based  on 
the  results  presented  in  References  17  and  18. 


- 6 - 


2.1.1,  Bfi8ic  Concepts 


The  lifting  surfaces  of  the  configuration,  (such  as  wing,  canard, 
horizontal  tail,  etc.) , are  divided  into  sub-surfaces  (Fig.  1).  Each 
subsurface  is  then  divided  into  trapezoidal  or  triangular  cells.  A 
bound  '*horse-shoe"  vortex  is  placed  in  each  cell,  as  described  in  Fig. 
1.  In  order  to  satisfy  the  boundary  condition  for  the  flow  to  be 
parallel  to  the  surface  of  the  wing,  a control  point  is  placed  in  each 
cell,  (Fig.  1).  The  shedding  of  the  free  vortices  (which  represent 
the  wake)  is  presented  as  an  option  for  the  user,  who  has  a large 
variety  of  possibilities  available  (for  examples  from  trailing  edge, 
side  edges,  leading  edge).  Each  free  vortex  is  divided  into  a number 
of  finite  segments,  and  a half-infinite  last  segment.  The  Biot-Savart 
law  enables  the  calculation  of  the  velocity  induced  by  a vortex  segment 
AB  of  strength  K on  a point  C in  the  flow  field: 


K • (cos  6ac  + cos  cTbX) 


(The  last  equation  is  in  scalar  form,  and  "r"  is  the  distance  between 
the  point  C and  the  segment  AB) . 

The  iterative  procedure  starts  with  the  calculation  of  a geometric 
influence  matrix,  with  the  free  vortices  being  straight  lines  shed  at 
an  angle  to  the  surface,  as  an  Initial  condition. 

The  boundary  condition  at  the  control  points  enables  the  construc- 


tion of  a set  of  linear  equations.  A typical  form: 
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Each  sub-matrix  represents  the  influence  of  subsurface  i 

on  subsurface  j.  (In  the  example  above  there  are  three  subsurfaces). 

The  order  of  each  submatrix  is  written  in  parenthesis,  and  n,  m,  k 
defines  the  number  of  cells  in  each  of  the  subsurfaces  1,  2,3,  according- 
ly. The  strength  of  the  bound  vortex  in  cell  number  j,  of  subsurface 
i is  denoted  by  K ^ , while  w*  is  the  component  of  the  free  stream 
velocity  normal  to  the  surface  of  the  wing,  at  the  control  point  of  cell 
nunfcer  j,  of  subsurface  i.  Column  vector  K*  is  to  be  solved.  When 
the  program  is  solved  for  the  strength  of  the  bound  vortices,  it  also 
calculates  the  aerodynamic  coefficients,  checking  for  their  convergence 
(except  for  the  first  cycle) . 

The  aerodynamic  coefficients  are  Cp,  Cp  , CM*  The  program  also 
calculates  the  center  of  pressure  X , according  to  the  following 
formulae : 
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The  lift  coefficient  (including  all  possible  terms) : 


N N 
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Here  AX^  , Ay^  , Az^  are  the  cartezian  components  of  a bound  vortex 
segment  ij  in  the  direction  of  x,  y,  z accordingly,  while  u^ , v^, 
w^  are  the  disturbance  velocities  calculated  in  the  middle  point  of 
each  vortex  segment  belong  to  the  wing  surface.  The  second  term  on  the 
right  hand  side  of  Eq.  (3)  is  of  second  order,  and  it  expresses  in  general 
form  the  part  of  the  in-surface  trailing  vortices  to  the  lift,  being 


usually  ignored. 


The  induced  drag  coefficient  and  the  pitching  moment  coefficient  are 


given  in  Eq.  (4)  and  (5): 
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, ZRe^  are  x and  z coordinates  of  pitching  moment  reference 
point.  The  second  term  on  the  right-hand  side  of  Eq.  (5)  expresses  the 
rolling  moment  because  of  induced  drag,  and  it  is  not  negligible  in  the 
calculations  of  complex  configurations,  such  as  canard-wing  at  high  angles 
of  attack. 


The  center  of  pressure  is  given  by  definition: 


X 

c.p. 


(6) 


(measured  from  X^^  ) 

In  general,  it  should  be  noted  that  some  of  the  assumptions  made 
when  neglecting  terms  like  those  containing  the  effect  of  sidewash, 
proved  to  be  incorrect  in  strong  interactions  between  canard  wake  and 
wing  surface,  and  in  the  case  of  strong  leading  edge  vortices.  Thus  the 
calculation  of  aerodynamic  coefficients  in  general  cases  should  include 
all  possible  contributing  factors. 


The  program  uses  an  iterative  process  for  the  re-calculation  of  the 
trajectories  of  the  free  vortices,  until  convergence  of  aerodynamic  co- 
efficients is  achieved.  (The  strength  of  those  is  calculated  after  the 
solution  for  the  strength  of  the  bound  vortices) . 


The  procedure  is  as  follows : 


Each  point  (i)  of  a segment  of  a free  vortex  is  moved  according  to 
Euler's  method  of  integration 
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Here  u^  , v^  *w^  are  the  Cartesian  components  of  the  induced  velocity 

(2)  (2) 

on  point  (i)  , before  moving  it.  , z ^ are  the  new  y and  z co- 

ordinates of  point  (i) . (The  flow  is  considered  steady). 

More  sophisticated  integration  methods  should  be  applied  in  cases  where 
strong  disturbances  are  expected,  such  as  in  the  cases  of  very  close  inter- 
actions between  lifting  surfaces,  or  very  low  aspect  ratios.  Hie  integration 
process  continues  iteratively  until  convergence  of  all  free  vortices  is 
achieved,  or  until  the  process  exceeds  a prefixed  number  of  iterations.  The 
influence  matrix  is  reclaculated , leading  to  a new  strength  of  each  bound 
vortex.  Aerodynamic  coefficients  are  then  recalculated,  and  all  the  process 
starts  again  until  convergence  of  the  coefficients  is  achieved  (see  flow- 
chart of  the  method,  Fig.  23). 

The  effect  of  compressibility  is  taken  into  account  by  applying 
Goethert  rule  and  incorporating  it  into  the  computer  program. 

The  accuracy  of  this  rule  and  its  application  to  the  vortex  flow  over 
wings  has  been  discussed  also  by  Rubbert  and  Brune  [30], 

By  using  this  compressibility  correction,  the  present  method  can  now 
be  used  over  the  range  of  subsonic  Mach  numbers  up  to  the  critical  Mach 
nunfcer . 

2.1.2.  Results  of  Calculations 
( a)  Conputational  Aspects 

The  use  of  the  singular  vortex  element  in  vortex  lattice  methods 
presents  some  difficulties.  The  main  one  is  the  too  higji  velocities  in- 
duced by  the  vortex  on  points  which  are  too  close  to  the  vortex  segment. 
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The  reason  for  that  Is  the  fact  that  a potential  singular  vortex  does  not 
exist,  due  to  viscous  effects.  Some  attempts  were  made  to  overcome  this 
problem,  for  example  Rom  and  Zorea's  [13]  attempt  to  use  the  model  of 
Sprelter  and  Sacks  [31]  for  the  calculation  of  the  radius  of  each  vortex 
core,  or  Maskew's  "Sub  Vortex  Techniques"  [1],[36],  One  of  the  main 
conclusions  from  these  attempts  was  that  It  is  satisfactory  In  moat  cases 
to  use  a general  "cut  off  distance"  for  all  the  vortices  In  the  model. 

The  last  method  saves  a lot  of  confuting  time,  and  it  was  used  In  the 
present  method.  The  results  show  that  it  is  reasonable  to  use  a "cut  off 
distance"  similar  to  the  width  of  the  partition  into  elemental  panels. 
Strong  vortices  such  as  those  shed  from  a leading  edge  may  require  a 
larger  "cut  off  distance". 

Another  problem  arises  when  shedding  too  many  free  vortices.  Because 
of  the  singularity  of  each  potential  vortex,  the  convergence  of  the  free 
vortices  trajectories  is  either  very  slow  or  completely  impossible.  The 
multiplicity  of  free  vortices  is  a result  of  a dense  partition  into 
elemental  panels,  (which  is  favorable  when  the  aim  is  to  achieve  a smoother 
load  distribution) . The  present  method  overcomes  this  difficulty  by  giving 
the  option  to  merge  close  free  vortices  into  one  vortex,  of  a strength 
which  is  the  sum  of  the  strength  of  the  merged  vortices.  The  nunfcer  of 
free  vortices  is  thus  reduced,  while  the  nuirfcer  of  elemental  panels  remains 
as  it  was . 


i 


Results  of  the  present  method  show  that  merging  too  close  free 
vortices  is  especially  effective  in  the  cases  of  leading  edge  separation, 
in  wings  of  low  aspect  ratio  (AR  < 2)  with  sharp  leading  edges.  In 


these  cases  the  convergence  was  much  more  rapid  and  the  nonlinear  aero- 
dynamic coefficients  calculated  showed  very  good  agreement  with  experimental 
results.  When  no  merging  was  applied,  the  calculated  results  showed  weaker 
nonlinear  behavior  than  expected,  resemling  experimental  results  for  wing 
with  round  leading  edges,  (see  for  example  Ref.  35).  The  effect  of  merging 
free  vortice  on  aerodynamic  coefficients  is  strong  when  dealing  with 
vortices  which  pass  upon  the  wing  surface;  otherwise  the  effect  is  neg- 
ligible. 

b)  Comparison  with  Experimental  Results 

1 . Closely  Coupled  Canard  and  Wing,  with  Leading  Edge  Separation 

The  configuration  was  chosen  for  comparison  with  Ref.  (32)  and  it  is 

shown  in  Fig.  2.  Both  the  canard  and  the  main  wing  are  of  aspect-ratio 

1.66,  and  of  similar  form.  Two  vertical  positions  of  the  canard  were 

checked,  as  it  is  shown  in  Fig.  2.  The  canards  are  vertically  located  at 

z/C  =0.2,  (-0.29)  above  (beneath)  the  main  wing  plane  accordingly,  and 
c 

their  aft  point  is  located  longitudinally  at  0.1Cc  from  the  main  wing 
vertex.  (Here  Cc  is  the  canard  mean  geometric  chord).  Behrbohm  [32] 
states  that  the  model  used  for  experiments  had  round  leading  edges,  which 
considerably  reduce  the  nonlinear  effects.  Indeed,  when  using  the  present 
method  with  the  option  of  merging  free  vortices,  (the  effect  of  which  was 
described  before  in  part  (a)  of  this  Chanter),  the  calculated  results  for 
the  lift  coefficient  of  the  high  canard  were  higher  than  the  experimental 
ones  (Fig.  3a,  upper  curve) . 

By  assuming  separation  of  vortices  without  merging  to  be  a better 
approximation  of  the  actual  separation  from  the  rounded  edges,  the  cal- 
culated results  obtained  were  better  than  those  previously  described, 


although  somewhat  lower  than  the  experimental  ones  (^Ig.  3a) . The  first 
results  however,  may  be  an  indication  for  the  maximum  nonlinear  effect 
that  can  be  obtained,  should  the  leading  edges  be  sharp. 

It  should  be  noted  that  the  present  method  also  gives  the  correct 
tendencies  for  the  low  canard,  in  comparison  with  the  high  canard.  The 
interesting  behavior  of  the  low  canard  can  be  observed  in  Fig.  3b,  where 
the  lift  coefficient  increases  linearly  up  to  an  angle  of  attack  of  approxi- 
mately 15°,  and  then  a strong  non-linear  effect  occurs  , tip  to  20°. 

This  phenomenon  has  been  observed  both  in  the  experimental  and  the 
calculated  results.  It  is  explained  by  the  fact  that  the  free  vortices 
of  the  canard  stay  below  the  main  wing  surface  up  to  a certain  point, 
and  as  the  angle  of  attack  increases  above  the  point,  the  free  vortices 
shed  from  the  canard  pass  above  the  main  wing  surface,  causing  positive 
interference  and  as  a result  non-linear  lift  is  obtained.  The  low  canard 
lift  coefficient  even  equals  the  high  canard  lift  coefficient  at  an  angle 
of  attack  of  20°. 

Correct  tendencies  for  the  high  and  low  canards  can  be  observed  also 
from  Fig.  4 where  the  change  in  the  pitching  moment  coefficients  is  descri- 
bed (versus  the  lift  coefficient). 

The  effect  of  interaction  between  high  canard  leading  edge  vortices 
and  wing  leading  edge  vortices  can  be  seen  from  Fig.  5.  Noteworthy  too 
is  the  fact  that  due  to  the  upwash  effect  of  the  wing  on  the  canard,  the 
leading  edge  vortices  of  the  canard  leave  its  surface  at  an  angle  higher 
than  half  of  the  angle  of  attack.  In  general,  this  accepted  initial  value 
which  was  fair  enough  for  single  wings,  mav  lead  to  numerical  problems  in 
some  of  the  strong  interaction  problems  such  as  the  case  of  a low  canard. 
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In  the  latter  case  numerical  problems  arise  because  of  unfitting  initial 
conditions,  when  the  wake  of  the  canard  tries  to  cut  through  the  wing 
surface.  The  low  canard  case  also  shows  tendencies  of  instability  in  the 
flew,  within  a limited  range  of  angles  of  attack.  However,  the  high 
canard  case  which  is  much  more  common  does  not  present  any  problems. 

An  example  of  the  load  distribution  is  presented  in  Fig.  6. 

2 . Canard  and  Wing  Interference  with  no  Leading  Edge  Separation 

The  configurations  chosen  were  three  variations  of  a basic  model 
experimented  by  Gloss  [33],  The  basic  form  of  the  configuration  is  shown 
in  Fig.  7.  The  model  consists  of  a trapezoidal  main  wing  of  AR  - 2.5 
and  a swept  back  canard  of  AR  = 4.12.  The  first  variation  (denoted  as 
configuration  No.  2)  consists  of  a 51.7  deg.  swept  back  canard,  which 
has  0 deg.  dihedral  angle.  The  second  variation  (denoted  as  confi- 
guration No.  3)  consists  of  a 51.7  deg.  swept  back  canard  which  has  18.6 
deg.  dihedral  angle.  The  third  variation  (denoted  as  configuration  No.  4) 
consists  of  a 60  deg.  swept  back  canard  which  has  0 deg.  dihedral  angle. 

All  canards'  roots  were  located  at  a height  of  z/C  = 0.185  above 
the  wing  plane,  (C  is  the  wing  mean  geometric  chord) , and  at  a longi- 
tudinal distance  of  x/C  = 1.304  from  the  model  nose  to  the  canard  lead- 
ing edge.  The  partition  chosen  for  the  calculations  was  NC  * 4 , NS  - 10 
(using  trapezoidal  cells)  both  for  the  canard  and  for  the  main  wing,  in 
all  cases. 

Free  vortices  were  shed  from  the  trailing  edges  and  from  the  side 
edges.  All  three  cases  showed  very  good  agreement  between  the  experimental 
and  the  calculated  results  as  can  be  seen  from  Figs.  8 to  12.  The  calcu- 


lated results  depart  from  the  experimental  results  only  when  stall  effects 


As  an  exangsle  of  this  phenomenon  a view  of  main  and  secondary  lead- 
ing edge  vortices  was  sketched  for  a delta  wing  of  AR  **  2,  at  20  deg. 
angle  of  attack  (Fig.  15).  The  calculated  relative  strength  of  those 
vortices  appear  in  a Table  (Fig.  15) . (Trailing  edge  vortices  were 
omitted  from  Fig.  15  in  order  to  get  a clearer  view)  . In  delta  wings 
of  AR  < 2 , the  appearance  of  the  secondary  vortex  is  found  to  be  at  an 
inner  region  than  the  example  given  as  indicated  by  experiments. 

4.  The  F-4-E  Planform  in  Compressible  Flow 


In  the  following  example  application  of  the  method  has  been  made  to 
analyse  a complex  planform  also  making  use  of  the  Goethert  rule  for  the 
evaluation  of  the  compressibility  effects. 

The  geometry  of  the  model  includes  the  main  wing  which  has  a dihedral 
of  12  deg.  on  its  outer  span  and  a "saw  tooth"  form  on  its  leading  edge. 
The  tail  planform  has  a negative  angle  fo  dihedral  (-23.25  deg.)  and  la 


operating  as  a mono-block  control.  Both  the  main  wing  and  the  tail  are 


start.  Trying  a denser  partition  of  NC  * 4 , NS  » 10  for  the  canard, 
with  NC  = 10,  NS  « 20  for  the  main  wing  had  very  little  effect  on  the 
results,  apart  from  the  higher  computer  time  needed  for  the  solution. 

An  interesting  comparison  is  shown  in  Fig.  13  and  14.  The  present 
method  is  compared  to  calculated  results  of  modified  Polhamus  methods 
(from  Rdf.  33).  It  can  be  seen  clearly  that  the  present  results  for  the 
lift  coefficient  are  in  better  comparison  with  the  experimental  results, 
in  addition  to  evaluation  of  the  load  distribution  and  the  strength  and 
form  of  the  wake,  which  are  not  calculated  by  the  modified  Polhamus 
methods . 

Vortex  Formation 


extended  to  the  center  line.  A schematic  description  of  the  geometry  of 
the  model  partitioned  into  elemental  panels  is  shown  in  Fig.  16.  The 
calculated  results  are  compared  with  the  experimental  results  of  Ref.  14. 

Very  accurate  results  are  obtained  for  the  lift  coefficient  variation 
with  the  angle  of  attack,  at  M = 0.6  and  also  at  M = 0.8,  up  to  the  stall 
conditions,  (Fig.  17). 

Good  results  are  also  obtained  for  tail  deflections  up  to  -16  deg. 

There  was  no  point  in  comparing  results  for  the  pitching  moment,  because 
of  the  difficulty  to  correct  the  calculations  for  the  presence  of  the 
complex  body;  however,  deflections  of  the  tail  show  the  right  tendencies 
in  the  calculated  pitching  moment. 

The  experimental  results  for  induced  drag  coefficient  fall  between 
the  calculated  results  for  the  case  of  full  leading-edge  suction  drag 
and  those  calculated  for  zero  leading-edge  suction  drag. 

c)  Cruciform  Configuration  in  Subsonic  Flow 

The  non-linear  characteristics  of  a cruciform  combination  of  rectangu- 
lar wings  are  evaluated  by  the  present  method  (at  0.0  Mach  nuutoer) . The 
configuration  is  combined  of  cruci formed  rectangular  main  wings  and  tails, 
described  in  Fig.  19.  The  calculated  lift  and  pitching  moment  coefficients 
show  a strong  non-linear  variation  with  increasing  angle  of  attack  (Fig.  20) . 

The  existence  of  the  "turning  point"  in  the  curves  of  C.  , Cw,  C_  (Fig. 

L M 

20)  is  explained  by  examining  the  path  of  the  free  vortices  shed  by  the 
wing.  At  about  12°  those  vortices  pass  very  near  and  belcw  the  upper  tail, 
causing  a considerable  loss  of  lift  on  that  part.  Further  increase  of  the 
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angle  ot  attack  would  cause  those  free  vortices  to  make  their  way  above 
tail,  gaining  lift  again.  It  should  be  noted  that  there  is  some  loss  of 
stability  in  the  region  of  the  turning  point,  as  seen  from  the  pitching 
moment  curve  slope.  Experimental  data  on  cruciformed  wings  configurations 
indicate  similar  variation.  Results  are  sensitive  to  the  relative  location 
and  shape  of  the  tail. 

d)  Computer  Time 

Computing  time  for  a complex  configuration  is  at  present  rather  high 
(about  2000  sec  on  IBM  370/168  type).  However,  computation  time  may  be 
considerably  reduced,  should  the  program  be  optimized.  The  use  of  the 
iterative  scheme  is  time  saving  in  itself. 


— 1 

- 18  - 

2.2.  A Method  for  Calculations  of  the  Aerodynamic  Characteristics 

of  Slender  Bodies  with  Vortex  Separation 

The  U9e  of  vortex  elements  enables  the  proposition  of  a simple 
model  for  representing  the  separation  of  vortices  from  the  body  surface. 

Mainly  the  model  of  the  body  is  that  of  an  "Annular  Wing",  built 
of  a vortex  lattice  (see  for  example  Ref.  2,23),  with  the  possibility  to 
separate  vortices  from  the  surface.  The  separation  point  in  the  model 
is  located  in  the  function  between  two  adjacent  bound  vortices  or  in 
some  point  located  on  the  line  of  the  trailing  in-surface  vorticity. 

The  vortex  arrangement  on  a cylindrical  body  is  shown  in  Fig.  21, 
together  with  an  example  for  the  method  of  separation  between  two 
adjacent  elemental  panels. 

This  method  of  separation  enables  also  to  separate  vortices  from 
pre-chosen  regions  of  the  bodv  surface,  according  to  the  separation 
theory  to  be  chosen. 

References  27,  28,  29,  46, for  example,  discuss  the  problem  of 
determining  a criterion  for  the  separation  of  vortices  from  the  surface 
of  slender  bodies  at  high  angles  of  attack.  However  a thorough  survey 
of  the  proposed  criteria  for  separation  is  under  consideration,  aimed 
to  choose  a general  model  which  should  be  also  as  simple  as  possible 
for  the  purpose  of  incorporation  with  the  method  of  calculation.  Ref. 

29  for  example,  presents  an  interesting  method  for  calculating  aero- 
dynamic characteristics  of  bodies  having  vortex  separation  at  high 
angles  of  attack,  but  the  cost  in  C.P.U.  time  for  utilizing  this 
method  seems  to  make  it  unattractive  for  incorporation  with  a general 
method  for  complex  configurations. 
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Other  potential  elements  such  as  point  sources  or  sinks  might  be 
needed  for  the  description  of  a blunt  nose  of  a body  for  example. 

Those  elements  are  distributed  inside  the  body  as  described  in  Fig. 

21. 

It  should  be  noted  that  the  calculation  of  wing-bodv  configuration 
is  carried  out  essentially  in  the  same  way  as  described  for  the  case  of 
wings  (Chapter  2.1.1.  in  this  report)  with  the  geometrical  influence 
coefficients  of  the  body  taking  proper  places  in  the  general  matrix 
(Eq.  2). 

A general  flow  chart  for  the  procedure  of  calculation  appears  in 
Fig.  23. 

A method  for  the  representation  of  thick  wings  appears  in  Fig.  22. 
Here  too  potential  elements  such  as  line  sources/sinks  might  be  needed 
for  a complete  description  of  thickness  and  of  round  edges,  (apart  from 
the  vortex  element) . 

For  more  accurate  calculation  of  the  local  pressure  coefficient^ 
methods  such  as  the  "Sub-Vortex  Technique"  by  Maskew  [1],[36]  may  be  used, 
after  adaption  to  the  three  dimensional  case. 

I 
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3.  EXPERIMENTAL  STUDY  OF  THE  AERODYNAMIC  CHARACTERISTICS  OF  SLENDER 
BODIES  AND  TOE  EFFECT  OP  INJECTED  JETS , AT  HIGH  ANGLES  OF  ATTACK . 

3.1.  Lateral  Forces  and  Moments  Acting  on  a Body  at  Hlgft  Angles  of 
Attack  at  Subsonic  and  Tranaonlc  Speeds 

The  occurrence  of  asymmetric  vortex  flow  and  hence  side  forces  and 
yawing  moments  on  bodies  of  revolution  at  high  angles-of-attack  and  zero 
sideslip  has  been  observed  experimentally  and  is  under  intensive  stddies 
at  various  laboratories . 

All  previous  investigations  have  shown  that  the  nose  shape  is  one 
of  the  most  significant  parameters  affecting  side  force  characteristics. 
It  has  been  found  experimentally  that  it  is  possible  to  reduce  side 
forces  by  the  use  of  devices  causing  boundary  layer  transition  or  vortex 
generators  on  the  nose. 

In  the  present  investigation  air  lets  blown  from  a body  of  revolu- 
tion are  used  in  order  to  reduce  the  side  forces. 

There  are  many  parameters  that  may  affect  these  side  forces  such 
as  the  position  of  the  stations  on  the  body  for  these  lets,  the  polar 
angle  of  each  station,  the  mass  flow  and  velocities  of  the  jets.  Some 
of  these  parameters  have  been  investigated  and  reported  in  this  paper. 

At  moderate  angles  of  attack  the  adverse  cross-flow  pressure  gra- 
dient on  the  leeside  of  the  bodv  causes  the  boundary  layer  to  separate 
and  then  to  roll-up  into  a system  of  distinct  vortices. 

The  rolled-up  vortex  sheet  may  stay  joined  to  the  body  and  be  con- 
tinuously fed  from  the  separated  boundary  laver  or  it  may  leave  the  body 
entirely  further  downstream.  Since  the  boundary  layer  separation  is 
caused  by  the  adverse  pressure  gradient,  the  separation  line  may  be 
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expected  to  be  close  to  the  line  of  the  minimum  pressure  coefficient  on 
the  leeward  side  of  the  body  [37]. 

Due  to  geometrical  irregularities  of  the  nose  of  the  body  or  irregu- 
larities in  the  flew  direction  one  side  of  the  boundary  layer  may  separate 
first  from  the  body,  retaining  a certain  vortex  strength  in  the  corresponding 
rolled-up  vortex  sheet,  the  other  side  of  the  boundary  layer  may  remain 
attached  and  separate  only  further  downstream  on  the  body  with  a corresponding- 
ly stronger  vortex.  Thus  we  find  an  asymmetric  vortex  system  generated  be- 
hind the  body  at  high  angles  of  attack.  This  is  supported  by  the  results  of 
Ref.  38  where  it  is  shown  that  the  side  forces  are  connected  to  asymmetry 
of  circumferential  pressure  distribution  and  asymmetry  of  vortex  sheet 
separation  on  the  leeslde  of  the  body.  In  the  visualization  tests  of  Ref. 

37  it  is  shown  that  as  the  angle  of  incidence  increases  the  separation  line 
is  shifted  towards  the  windward  side.  Hie  circumferential  angles  of  sepa- 
ration ^3i*^82^  as  ™*asured  from  the  windward  meridian,  indicated  a 
non-symmetrical  separation  on  the  cylindrical  afterbody  when  laminar  boundary 
layer  conditions  exist.  Hiis  asymmetry  increased  with  increasing  angle  of 
attack.  The  maximum  angular  difference  very  nearly  coincided  with  the 
maximum  measured  side  force.  (The  maximum  angular  difference  was  of  the 
order  of  20°  to  40°  for  the  sharp  and  slightly  blunted  slender  noses 
respectively,  while  this  difference  decreased  to  5°  to  10°  for  the  lowest 
fineness  ratio  nose) . 

At  high  angles  of  attack  the  side  vortices  could  grow  and  separate 
at  several  points  along  the  body:  this  phenomenon  was  found  by  [37], 

[38]  and  [39]. 
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The  phenomenon  of  these  lateral  forces  can  be  also  modeled  by  the 
urapulsively  started  cross-flow  analogy.  The  calculation  of  side  force* 
and  pressure  distributions  on  Inclined  cylindrical  bodies  at  high 
angles  of  attack  is  presented  by  Lamont  and  Hunt  in  Ref.  40  and  41. 
Using  the  cross-flow  analogy  the  induced  lateral  flow  field  is  evaluted 
without  specific  reference  to  the  state  of  the  boundary  layer  on  the 
body.  The  validity  of  this  assumption  must  be  studied. 

The  generally  accepted  results  of  the  investigations  of  the  pheno- 
mena associated  with  vortex  asvmmetry  on  slender  bodies  at  high  angles 
of  attack  are  the  following: 

(1)  The  initial  direction  of  the  side  force  i9  unpredictable  because 
of  its  connection  with  small  irregularities  in  the  nose  geometry; 
however  once  the  direction  is  established  this  does  not  change  as 
the  incidence  angle  is  increased. 

(2)  Reynolds  number  has  no  effect  on  the  angle  of  onset  of  asywnetry 

but  does  have  an  effect  or  the  magnitude  of  the  side  force. 

(The  Reynolds  number  effect  is  in  defining  the  boundary  layer  in 

the  sense  of  whether  it  i9  supercritical  or  subcritical.  For 

Instance  in  Ref.  40  the  relation  of  in  the  supercritical  to 

max 

the  subi.tii.ical  is  1.5/2. 5). 

(3)  The  magnitude  of  the  side  forces  increases  with  increasing  the 
finess  ratio  of  the  nose. 

(4)  Side  forces  remain  small  below  a - 25°.  Beyond  this  they  in- 
crease and  reach  peaks  st  around  35°  - 40°,  depending  on  the  Mach 
number  and  geometrical  configuration. 


(5)  The  magnitude  of  the  side  force  is  decreased  by  as  much  as  80% 

when  a boundary  layer  strip  is  installed  on  the  windward  side. 

This  is  a result  of  early  transition  to  turbulent  boundary  layer. 

Side  forces  and  moments  may  be  potentially  hazardous  to  the  control 
and  stability  of  slender  configurations  such  as  fighters  and  modem 
missiles,  at  high  angles  of  attack.  These  side  forces  and  moments  may 
be  overcome  by  sufficient  control  authoritv  and  autopilot  gain  margin 
so  that  flight  control  is  maintained,  or  by  aerodynamic  devices  to 
suppress  the  asymmetric  vortex  pattern. 

In  Refs.  [42,43  and  44]  the  investigators  have  pointed  out  that 
the  nose  shape  is  a most  significant  parameter  affecting  side  forces, 
therefore  adding  transition  strips  or  vortex  generators  on  the  nose 
can  reduce  the  side  forces . 

These  devices  cause  regular  svmmetric  boundary  laver  flow  on  the 
nose  of  the  body,  thus  eliminating  the  previously  discussed  asymnetric 
boundary-layer  separation  and  so  the  associated  vortex  system  is  sym- 
metric, therefore,  eliminating  the  large  side  forces  and  moments. 

In  the  work  of  Keener  and  Chapman  [42]  reduction  of  side  forces  is 
achieved  by  adding  meridianal  strips  of  grit  at  the  angular  position 
A ■ i 30°.  In  Ref.  [43],  the  positioning  of  a ring  on  the  nose  caused  a 
great  reduction  in  the  large  vawing  moments  at  high  angles  of  attack. 

In  the  work  of  Clark,  Peoples  and  Briggs  [44]  the  additon  of  vortex 
generators  at  the  nose  was  sufficient  to  cause  the  separated  vortices 
to  be  more  symmetric  at  the  critical  angles  of  incidence.  In  this 
present  investigation  air  lets  blown  from  the  bodv  of  revolution  are 
used  to  affect  flow  svnnetry  and  reduce  the  side  forces  and  moments. 
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This  is  shown  by  an  experimental  investigation  wherein  the  side  forces 
and  moments  are  measured  with  the  air  lets  blown  at  different  blowing 

coefficients,  C . 

u 

Preliminary  tests  have  been  carried  out  in  the  high  subsonic  and 
the  transonic  flow  regions . 

3.2.  The  Effects  of  Jets  Injected  Near  the  Nose  of  the  Body  on  the 

Lateral  Forces  and  Moments  at  High  Angles  of  Attack,  in  Subsonic 
and  Transonic  Speeds. 

3.2.1.  The  Model  and  Facilities  for  the  Experiments 

The  experiments  described  herein  were  conducted  in  the  Subsonic 
Wind  Tunnel  of  the  Aeronautical  Research  Center  of  the  Technion  with  a 
cross  section  of  lm  x lm. 

The  model  shown  in  Fig.  24  is  an  axisvmmetric  body  of  revolution 
with  a body  diameter  of  3 in.  and  an  overall  finess  ratio  of  6. 

Three  different  forebodies  could  be  used;  a pointed  cone,  a blunt 
cone  and  a blunt  ogive  (shown  in  Fig.  24) . 

After  conducting  some  preliminary  tests  it  was  decided  to  blow  lets 
from  two  sections,  one  very  close  to  the  nose  apex  and  the  second  further 
along  the  nose  (Fig.  24) . At  each  station  three  pairs  of  holes  were 
drilled  at  ± 60°  and  ± 30°  as  indicated  in  Fig.  24  and  the  diameter 
of  the  inlection  holes  was  1.2  run  perpendicular  to  the  bodv's  axis 
(this  diameter  was  chosen  so  as  to  obtain  blowing  velocities  of  the 
same  order  of  magnitude  as  the  tunnel  flow  speed)  . 

A system  of  rigid  and  flexible  tubing  was  arranged  in  the  model  so  as 
to  obtain  svmmetrical  blowing  of  the  lets  through  a different  pair  of 
holes  each  time. 
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Throughout  the  experimental  program,  a six  component  internal  strain 
gage  balance  was  used  to  measure  the  aerodynamic  forces  and  moments.  This 
balance  was  very  sensitive  to  side  forces. 

The  experimental  system  is  capable  of  changing  the  rate  of  flow  and 
of  providing  measurements  of  the  pressures  and  of  the  rates  of  flow. 

The  experiments  were  conduced  at  a flow  speed  of  20m/sec  (for  a 


Reynolds  number  of  1.08  x 10  based  on  the  diameter)  with  the  pointed  cone 
configuration  and  for  a range  of  angles  of  attack  of  30°  - 56°.  In 
the  preliminary  tests  we  found  this  configuration  to  give  the  highest 
values  of  side  force. 

Tests  with  no  injection  and  a few  tests  with  injection  were  re- 
peated in  order  to  verify  the  repeatability  of  the  test  results.  The 
jets  were  also  blown  without  flow  and  the  measured  side  force  was  found 
to  be  negligible,  showing  that  the  blowing  is  symmetrical. 

Qualitative  tests  for  flow  visualization  were  also  performed. 

Using  oil  paints  on  the  model,  the  development  of  the  stream  lines  on 
the  body  was  followed.  The  flow  patterns  indicated  positions  of  sepa- 
ration as  shown  in  Fig.  25. 

3.2.2,  Results 

The  results  of  the  tests  can  be  summarized  according  to  the  angular 
position  of  the  lets  on  the  model: 

Blowing  at  60°-  At  both  stations  increased  values  of  the  side  force  co- 
efficients in  comparison  to  those  obtained  without  blowing  were  observed. 

Blowing  at  + 30°  - At  station  II,  high  values  of  side  force  coefficients 
were  observed  with  increasing  rate  of  flow.  At  station  I,  two  interesting 
results  were  observed  (see  Fig.  26)  - the  line  of  C - 166.33  x 10-\ 
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where  the  jets  kept  the  side  force  coefficient  from  growing  until  40° 

and  then  caused  a peak  higher  than  the  original;  and  the  experiments  for 

C » 289.34  x 10  ^ , in  which  the  side  force  coefficient  followed  the 

u 

original  line  to  an  angle  of  40°  where  there  was  a change  of  sign  until 
a negative  peak  was  reached  at  44°,  followed  by  a rise  to  a peak  similar 
to  the  original. 

Blowing  at  30  At  station  II,  a reduction  of  the  maximum  side  force 

coefficient  up  to  50°  was  observed  in  the  two  experiments  shown  in 
Fig.  27. 

Both  the  curves  C = 33  25  x 10  and  285.'/  x 10  ^ d :crcas*-: 

l 

from  their  maxima  at  a = 46°  to  zeio  side  force  coefficient  at  a * 53°. 

The  yawing  moment  coefficient  was  decreased  in  these  experiments 
by  up  to  30%, 

At  station  I,  a meaningful  reduction  of  side  force  for  certain  ang^r.- 
of  attack  was  observed  (Fig.  28) . 

On  the  yawing  moment  coeffcient  curve  (Fig.  29)  we  observed  a meaning- 
ful reduction  with  jets  blowing  compared  to  the  original  curve 

On  examination  of  the  side  force  coefficient  versus  the  rate  of  flow 
diagram  (Fig.  30)  it  is  obvious  that  there  is  an  optimum  rate  of  flow  for 
blowing  the  lets  . The  arrows  on  the  Figure  point  to  a concentration  of 
many  test  results  (on  this  curve  the  side  force  direction  is  opposite 
to  its  direction  on  Fig  31) . 

Side  Forces  Versus  Reynolds  Number 

The  experiments  of  this  investigation  were  conducted  in  a iange  of 
Reynolds  number  from  1.07  x 10^  to  1.69  x 10^ , based  on  the  body’s 


diameter . 
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Fig.  31  shows  the  change  of  side  force  coefficient  versus  crossflow 
Reynolds  number,  at  = 21.6  x 10  ^ . 

Consistency  of  Results 

Fig.  32,  which  describes  the  experiments  without  blowing  shows  good 
consistency  up  to  48°.  From  this  point  scattering  of  results  can  be 
observed  but  they  are  still  acceptable  as  consistent  (the  scattering  of 
results  is  a characteristic  phenomenon  in  all  works  on  this  subject. 

It  is  common  practice  to  attribute  this  to  the  fact  that  the  side  forces 
are  sensitive  to  changes  in  the  nose  alignment  and  in  the  experimental 
conditions  and  also  to  mounting  induced  vibrations. 

In  the  blowing  experiments  which  were  repeated,  similar  amplitudes 
were  observed  but  the  sign  of  the  side  force  changed  in  some  cases. 

From  the  experiment  that  was  conducted  without  flow  it  was  found  that 
no  side  forces  were  measured  by  the  balance  and  the  conclusion  is  that  the 
blowing  during  the  experiments  was  symmetrical. 

Analysis  of  Oil  Flow  Photographs 

From  a comparison  of  the  pictures  with  no  blowing  to  those  with 
blowing  which  was  carried  out  for  two  angles  of  attack,  a » 46°  and 
a = 48°  (on  which  the  maximum  no  blowing  side  force  coefficient  was 
observed) , the  following  conclusions  can  be  derived: 

(1)  The  blowing  of  the  jets  delays  the  separation  point  to  much 
further  downstream  on  the  body  (because  of  the  transition  to 
turbulent  flow)  . 

(2)  Hie  flow  separation  lines  on  each  side  become  more  symmetrical. 

These  results  are  similar  to  the  results  described  in  Refs .40 ,41 
and  42,  in  which  the  investigators  added  rigid  elements  on  the  nose. 
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Preliminary  Tests  in  the  High  Subsonic  and  in  the  Transonic  Regions 

The  teats  were  carried  out  in  the  transonic  wind  tunnel  of  the 
Aeronautical  Research  Center-  of  the  Technion,  which  has  a cross  section 
of  60cm  x 80cm.  The  model  is  the  pointed  cone  nose  (Fig.  24) . The 
investigations  included  experiments  with  and  without  let  injection  at 
Mach  numbers  of  0.6  and  0.8,  for  angles  of  attack  between  35°  to  56°, 

The  injection  took  place  at  station  A (Fig.  24) , at  an  angle  of  -30° 
and  with  a blowing  coefficient  of  - 149.4  x 10  7 for  M * 0,6 
and  q,  = 84.07  x 10-7  for  M = 0.8. 

Schlieren  photographs  were  taken  during  the  tests.  At  M * 0.6 
the  blowing  had  only  little  effect  on  the  reduction  of  the  side  force 
coefficient  (Fig.  33) , although  it  reduced  to  some  extent  the  yawing 
moment  (Fig.  34).  However,  at  M = 0.8  both  side  force  and  yawing  moment 
were  considerably  reduced  within  the  range  of  the  angles  of  attack 
35°  < a < 48°  (Figs.  35,  36).  The  effectiveness  of  the  blowing  disappeared 
for  angles  of  attack  higher  than  48°. 

From  Shlieren  photographs ( Fig . 37)  it  can  be  seen  that  a shock  wave 
appeared  near  the  shoulder  of  the  body  at  M = 0.8,  However,  the  flow 
picture  with  jet  blowing  (Fig.  37b)  is  different  from  the  flow  without 
blowing  (Fig.  37a) . The  considerable  reduction  in  side  force  and  yawing 
moment  at  M = 0.8  encourages  further  investigation  for  the  effect  of 
jet  blowing. 

Effects  of  Reynolds  number,  influence  of  the  rate  of  blowing  and 
shock  wave  formation  should  be  investigated  at  these  and  higher  transonic 
Mach  numbers . 
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4 . DISCUSSION  AND  CONCLUSIONS 
4.1 

The  present  method  offers  an  efficient  means  of  analysis  and  design 
of  complex  aerodynamic  configurations.  The  good  results  obtained  for 
strong  interaction  cases  and  for  non-linear  effects  encourage  further 
development  and  refinement  of  the  method.  One  aim  in  future  development 
has  been  mentioned  in  Section  2.2.  of  this  report,  and  that  is  including 
the  presence  of  bodies  with  vortex  separation.  The  presence  of  a body 
is  influencing  especially  the  pitching  moment  coefficient  of  wing-body 
configurations,  and  also  the  form  of  the  wake,  the  pressure  distribution 
and  the  total  lift  coefficient  at  high  angles  of  attack.  Another  aim 
in  future  development  should  be  inclusion  of  thickness  effects  for  lifting 
surfaces . 

Although  expensive  for  use  in  the  sense  of  computer  time  and  memory, 
a correct  description  of  thick  wings  with  round  leading  edges  might  give 
more  correct  pressure  distribution  and  might  enable  a good  estimate  of 
the  Induced  drag,  and  of  the  leading  edge  suction  force. 

Thought  should  be  given  also  to  get  "smoother"  pressure  distributions, 
although  present  results  show  that  the  total  aerodynamic  coefficients  are 
in  good  comparison  with  experimental  results  even  with  the  simple  "block" 
distribution  used. 

It  should  be  concluded  that  the  present  method  can  be  used  for  the 
calculation  of  aerodynamic  characteristics  of  complex  configurations  at 
high  angles  of  attack  and  at  subsonic  speed,  with  reasonable  cost  in 


The  Modified  Vortex  Lattice  Method  for  Calculations  of  Cot 
Wing  Planformls  at  High  Angles  of  Attack . 
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computing  time  due  to  the  simplicity  of  the  basic  model  used. 

Using  the  present  method,  the  analysis  holds  up  to  the  effects  of 
stall  or  vortex  break-down. 

4.2.  The  Effects  of  Jets  Injection  Ne.ar  the  Nose  ot  a Slender  Body  on 

the  Side  Kn and  Yawing  Moment  at  High  Angles  of  Attack 

From  the  zesults  of  this  research  it  can  be  seen  that  it  is 
possible  to  reduce  side  forces  and  even  change  their  sign  by  blowing 
symmetric  lets  from  the  body  as  close  as  possible  to  the  nose  apex 
(the  importance  is,  of  course,  in  reducing  the  maximum  side  forces). 

The  most  effective  results  appear  on  blowing  at  an  angle  of  - 30° 
(on  the  windward  side)  . It  can  be  seen  that  the  let  blowing  is  inef- 
ficient at  angles  of  attack  above  52°. 

The  amount  of  blowing  used  did  not  cause  any  change  in  the  normal 
force  drag. 

The  rate  of  flow  of  the  jets  needed  to  reduce  the  side  forces  is 
very  small.  This  solution  is  not  sensitive  to  changes  in  the  rate  of 
flow;  increasing  the  rate  of  flow  did  not  cause  significant  changes  in 
side  force  coefficient. 

Effect  of  Angle  of  Attack  - Increasing  angle  of  attack  causes  the  side 
force  to  increase  until  it  reaches  a peak  at  about  50°.  Above  this 
angle  further  increase  of  angle  of  attack  causes  decrease  of  the  side 
forces  . 

Angle  of  Onset  Asymmetry  - From  the  results  of  the  experiments  it  is 
clear  that  the  blowing  iets  do  not  cause  any  change  in  the  value  of  the 
angle  of  onset  of  asymmetry. 
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Jet  Control  - More  effective  results  in  utilizing  blowing  jets  may  be 
achieved  by  designing  a control  system  which  changes  the  rate  of  flow 
as  a function  of  angle  of  attack. 

A control  system  can  be  devised  to  vary  the  rate  of  flow  in  order 
to  keep  low  values  of  side  force  coefficients. 

Effect  of  Reynolds  Number  and  Mach  Number  Changes 

At  this  stage  of  tests  it  is  impossible  to  define  the  effects  of 
Mach  number  on  side  forces. 

l urther  tests  are  required  in  order  to  investigate  the  effects  of 
both  Mach  number  and  Reynolds  number  on  the  lateral  forces  and  moments 
of  these  slender  configurations  at  high  angles  of  attack  with  and 
without  blowing. 
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FIGURE  1 - DESCRIPTION  OF  THE  MODEL. 


PIGURE  2 - COMPARISON  WITH  EXPERIMENTAL  RESULTS  - REP. 

CANARD  CONFIGURATION  NO.  1 PARTITIONED  INTO 
ELEMENTAL  PANELS. 
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PAIRED  DATA  CURVE  (ROUND  L.E.),HIGH  CANARD 

PRESENT  CALCULATION 

MAX.  LIFT  PREDICTED  BY  PRESENT  WORK  FOR  SHARP  L.Z. 


LINEAR  THEORY 

(EXPERIMENTAL  DATA  TAKEN  FROM  REF. 32) 
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PIGURE  3-a  - LIFT  COEFFICIENT  Vs.  ANGLE  OF  ATTACK 
FOR  CONFIGURATION  No.  1,  HIGH  CAHARD. 


4 

4 


- 40 


PAIRED  DATA  CURVE  (ROUND  L.E. ) ,LQW  CANARD 

PAIRED  DATA  CURVE  (ROUND  L.E.),HIGH  CANARD 

PRESENT  CALCULATION,  LOW  CANARD 

LINEAR  THEORY 
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FIGURE  3-b  - LIFT  COEFFICIENT  Vs.  ANGLE  OF  ATTACK 
FOR  CONFIGURATION  No.  1,  LOW  CANARD. 
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FIGURE  7 - COMPARISON  WITH  EXPERIMENTAL  RESULTS  - REF.  33  , WING  CANARD 
CONFIGURATIONS  Nos.  2,  3,  4. 
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FIGURE  9 


LIFT  COEFFICIENT  Vs.  DRAG  COEFFICIENT,  CONFIGURATION  No.  2. 
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FIGURE  10  - PITCHING  MOMENT  COEFFICIENT  Vs.  ANGLE  OF  ATTACK,  CONFIGURATION  No.  2 
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FIGURE  14  - COMPARISON  OF  EXPERIMENTAL  RESULTS  TO  VARIOUS  THEORETICAL 
RESULTS  FOR  CONFIGURATION  No. 2. 
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FIGURE  16  - SCHEMATIC  DESCRIPTION  OF  THE  MODEL  OF  AN  F-4-E 
PARTITIONED  INTO  ELEMENTAL  PANELS. 
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FIGURE  17  - LIPT  COEFFICIENT  Vs.  ANGLE  OF  ATTACK,  F-4-E  MODEL, 
AT  0.6  MACH  NUMPER. 
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LIFT  AND  PITCHING  MOMENT  COEFFICIENTS  Vs.  ANGLE  OF  ATTACK;  INDUCED  DRAG 
COEFFICIENT  Vs.  LIFT  COEFFICIENT,  CALCULATED  FOR  THE  CRUCIFORM  CONFIGURATION 
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FI  CURE  23  - SCHEMATIC  FLOW-CHART  OF  THE  CALCULATION , (PRESENT  METHOD) 
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FIG.26-  SIDE  FORCE  COEFFICIENT  VERSUS  ANGLE 
OF  ATTACK  FOR  INIJECTION  AT  STA.  I ANGLE  .30* 

( POINTED  CONE  ) 
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FIG. 28-SIDE  FORCE  COEFFICIENT  VERSUS  ANGLE 
OF  ATTACK  FOR  INJECTION  AT  STA.  I ANGLE  - 30* 
(POINTED  CONE  V = 20  — ) 


FIG.29- YAWING  MOMENT  COEFFICIENT  VERSUS  ANGLE 
OF  ATTACK  FOR  INJECTION  AT  STA.  I ANGLE  -30* 
(POINTED  CONE,V=20  ) 


FIG.  30-  CHANGE  OF  SIDE  FIG. 31-SIDE  FORCE  COEFFICIENT 

FORCE  COEFFICIENT  AT  46*  ANGLE  VERSUS  CROSS  FLOW  REYNOLDS 
OF  ATTACK  FOR  INJECTION  AT  STA.  NUMBER  AT  a =45*  Cm  =21.6*  10‘! 
TANGLE  -30*.(  POINTED  CONE) 


FIG. 33 -SIDE  FORCE  COEFFICIENT  Vs.  ANGLE  OF  ATTACK, 
AT  0.6  MACH  NUMBER  Cp.  =149. 4 * 10'7 
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AT  0.6  MACH  NUMBER  Cu  =149.4  * 10'7 
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FIG.  37 -SCHLIEREN  PHOTOGRAPHS  OF  THE  BODY  AT  0.8 

MACH  NUMBER,  a =44° 


